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ABSTRACT 

We have examined the statistics of the polarization position angles deter- 
mined for point sources in the NRAO-VLA sky survey (NVSS) and find that 
there is a statistically significant bias toward angles which are multiples of 45 
degrees. The formal probability that the polarization angles are drawn from a 
uniform distribution is exponentially small. When the sample of those NVSS 
sources with polarizations detected with a signal to noise >3 is split either 
around the median polarized flux density or the median fractional polariza- 
tion, the effect appears to be stronger for the more highly polarized sources. 
Regions containing strong sources and regions at low galactic latitudes are 
not responsible for the non-uniform distribution of position angles. We iden- 
tify CLEAN bias as the probable cause of the dominant effect, coupled with 
small multiplicative and additive offsets on each of the Stokes parameters. Our 
findings have implications for the extraction of science, such as information 
concerning galactic magnetic fields, from large scale polarization surveys. 
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1 INTRODUCTION 

Low frequency measurements of the polarization properties of point sources have been used 
to deduce information concerning galactic and extragalactic magnetic fields with some suc- 
cess (see, for example, Carilli & Taylor, 2002; Beck, 2001, and references therein). Future 
large scale projects such as the Square Kilometre Array (SKA) intend to make measurements 
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of ~ 2 X 10^ compact polarized extragalactic sources (see discussion in Carilli & Rawlings, 
2004) in order to understand the nature of the galactic magnetic field structure (Beck & 
Gaensler, 2004) and various aspects of cosmic magnetism (Feretti, Burigana & Ensslin, 2004; 
Feretti & Johnston-HoUitt, 2004). 

The statistics of polarization position angles (sometimes used in conjunction with struc- 
tural/jet position angles for spatially resolved sources) have also been used to investigate the 
global rotation of the universe. In particular, there is a history of claimed detections of global 
anisotropy and rotation effects using these methods (Birch, 1982; Phinney & Webster, 1983; 
Kendall & Young 1984), the most recent being made by Hutsemekers et al (2005) on the 
basis of the optical measurements of quasars with magnitudes ~ 14 — 15. However, see Joshi 
et al (2007) for a discussion of these effects in the context of the radio polarizations measured 
for flat spectrum sources in the 8.4GHz JVAS/CLASS survey (Jackson et al, 2007). 

The veracity of these interesting applications requires the data to be free of systematics to 
an unprecedented level of accuracy. In any survey of ~ 2 x 10^ sources there is the possibility 
for systematic effects at the level of 0.02% to have an impact on some statistics computed 
from the data. The aim of this paper is to investigate the systematics in the measurement of 
polarization in the NRAO-VLA sky survey (NVSS), the largest survey of point sources with 
published polarization information presently available (Condon et al, 1998). To our surprise, 
for such a widely used survey, wc have identified an obvious and statistically significant bias 
in the measured polarization position angles. 

Details of the NVSS survey are given in Condon et al (1998). Here we briefly summarize 
some of the survey details which are relevant to the following discussion. NVSS is a blind 
continuum survey which covers 82% of the celestial sphere at a frequency olv — lA GHz in 
Stokes parameters /, Q and U . It was made with the Very Large Array (VLA) operating in its 
D configuration. A total of 2 x 10® discrete sources were found with fiux densities / > 2.5mJy 
and the average noise levels were cx/ ~ 0.45 mJy beam~^ and ctq = cr^/ ~ 0.3 mJy beam~^. 
The sky was covered using 217446 overlapping, snapshot images. Each snapshot image was 
CLEANed (Hogbom 1974) using super-uniform weighting to optimise the beamshape. In the 
process of CLEANing large fields, the effect of CLEAN bias (section 3.1) can in principle 
reduce the flux density in the flnal images. In the NVSS survey, the CLEAN procedure for 
the total-intensity maps was terminated when the peak residual in the residual map fell to 
0.75 mJy bcam~^. Sources were extracted from the total intensity maps using an elliptical 
Gaussian-fitting procedure. In the case of the polarization (Q and U) images, a small amount 
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of CLEANing was carried out, and the flux densities were extracted by interpolation at the 
fitted total-intensity positions rather than by direct fitting. 

In section 2 we identify and characterize what we believe is a previously unidentified sys- 
tematic in the polarization position angles in the pubhshed NVSS catalogue. In sections 3.1 
and 3.2 we discuss the effects of CLEAN bias and offsets on the polarization angle distribu- 
tion and in section 3.3 we show how a combination of these effects can be used to explain 
the observed biases in the NVSS data. We make an attempt to "correct" these biases using 
a simple prescription in section 3.4. In the final section we discuss the possible implications 
of our findings for future observations. 

2 STATISTICS OF POLARIZATION POSITION ANGLES IN NVSS 

We have selected a sample of 519713 sources from the NVSS catalogue with flux / > lOmJy 
for which there is a quoted polarization position angle a defined by 



which is in the range —90° < a < 90°. We have computed a histogram of these data with 
bin widths of 5° and 10° (36 and 18 bins respectively); the results are presented in Fig. 1 
(top). It is clear from this that there is a substantive systematic bias toward position angles 
which are multiples of 45° and which correspond to states with either pure Q or pure U. 
The amplitude of the bias is around 5-10% whereas the expected random error on each 
(~ iV-^bin) is less than 1% in both cases. Since the effect is coherent over a number of 
bins it is extremely significant. Using a simple f^st one can compute the probability that 
the observed values of a are selected from a uniform distribution. In both cases the formal 
probability is < 10"^'^. Clearly, the null hypothesis, a uniform distribution, can be excluded 
at extremely high significance. 

The polarized flux is not detected at high signal-to-noise {S/N) in a large fraction of the 
sources and, therefore, we have for much of the subsequent discussion excluded sources with 
polarized flux density P < ImJy ; the resulting sub-sample contains a total of 203097 sources. 
Given the noise levels quoted for the survey this corresponds to a sample with a detection 
of polarization with S/N > 3. The equivalent histogram for this sample is presented in 
Fig. 1 (bottom). The bias noted for all sources is still detected at high significance in this 
sub-sample. 




(1) 



4 R.A. Batty e et al 




2.95xl0< - 




1 



n 



r 



If 




Figure 1. Histogram of the polarization position angle for all NVSS sources with flux greater than 10m Jy and a quoted angle 
(top) and those with P > Imjy (bottom), using 36 bins (left) and 18 bins (right). Note the apparent bias toward angles which 
are multiples of 45° which is present at extremely high significance in each. 



Such a bias should be regarded with some suspicion since a number of possible systematic 
errors might cause such an effect. Therefore, we have made a number of cuts on the data: 

(i) Examination of Fig. 2 which presents the distribution of polarization angles on the sky for 
/ > 300mJy illustrates that there are regions, that appear to be connected with the positions 
of known strong sources (for example, the Cygnus region) and low galactic latitude, where 
the polarization position angles are strongly correlated. It is conceivable that some of the 
catalogued sources in regions around strong sources are not actual objects, but are created 
by the aliasing power from the nearby single strong source via the sidelobes. Moreover, a 
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Figure 2. The polarization position angles of all NVSS sources with I > 300mjy using a Hammer-AitoflF projection. Each 
source is represented by line at the relevant angle to the vertical. The low galactic latitude plane is clearly visible as is the 
Cygnus region (RA~ 20hr, Dec~ 40°). 

similar effect might also be present due to strong emission from the galactic plane/centre. 
We would like to be sure that these sources are not responsible for the bias presented in 
Fig. 1. In order to investigate this we have made two cuts, the first which should exclude 
all these regions and the other which is very conservative: (a) we have excluded sources 
which are within 4° from sources with flux density / > 8Jy and also those within 10° of 
the galactic plane {\b\ < 10) and (b) exclude regions within 4° of source with flux density 
/ > 2Jy and those with |6| < 30°. The results are almost identical to the case when all 
sources are included and we conclude that aliasing from strong sources is not responsible for 
the bias. 

(ii) We have split the sub-sample of sources with polarization S/N > 3 about the median 
polarized flux Pmed(= y/Q^~+TP) ~ 2.1mJy and about the median polarized fraction nined(= 
P/I) ~ 0.066 and the results are presented in Fig. 3. In this flgure we have used the fractional 
differences N — {N — N)/N from a constant mean N are used instead of the total number 
N, since the samples are of different size. The level of fluctuations in the distribution for 
P > Pmed is a little higher than for the whole sample, whereas for P < Pmed it is more 
uniform, but still shows some weak evidence for the effect. This is to be expected since the 
low P sample will be affected by noise, which will have the tendency to make the distribution 
of a more uniform, even in the presence of a systematic effect. The histogram for 11 > limed 
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Figure 3. Splits of the data illustrating the apparent independence of the effect presented in Fig. 1. (left) the sample split 
about the median polarized flux P^cd ~ 2.1mjy; (right) the sample split about the median polarization fraction Hmcd = 0.066. 
In both cases the top is the sample above the median and the top that below the median. Note that we have used the fractional 
difference from the mean, TV, in these figures as opposed to the number in each bin N which was used in Fig. 1. 



is similar to that for P > Pmed? while that 11 < limed still appears to exhibit similar biases, 
albeit at a lower level. It is remarkable that the bias appears to be strongest in the regime 
P > Pmed and n > limed where the measured values of Q and U should be most reliable. 

On the basis of the discussion above we conclude that there are biases in the polarization 
position angles of sources which are in the NVSS catalogue. The bias toward multiples of 
45° suggests some kind of instrumental origin. Moreover, any effect which is present in the 
sky, such as those discussed in Hutsemekers et al (2005) and Joshi et al (2007), would have 
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been suppressed by Faraday rotation at this relatively very low frequency. In subsequent 
sections we will attempt to explain these data in terms of various systematic effects. 

Since the biases in the polarization angle histograms appear to be periodic, it seems 
interesting to characterize them using a phenomenological model created from sinusoidal 
waves with different periodicities. In sections 3.1 and 3.2 we will argue that there are three 
separate harmonics involved each of which is a result of a very different kind of systematic 
that one might expect to be present in the data at some level. In particular we will write 

the fractional difference from the mean as 
N - N 

N = — — ^ Acos8a + B cos Aa + C cos 2{a - ccq) , (2) 

where A, B, C and arc parameters. 

The main effect in the data has period 45° corresponding to the maxima at —90°, —45°, 
0°, 45° and 90° in Fig. 1 which is represented by the cos So; term in (2). In addition to 
this there appears to be a modulation of the peak heights such that those at —90°, 0° and 
90° are higher than those at —45° and 45°. This is represented by the cos4q; term. Both 
these harmonics seem to be present both in the whole sample and in the sample with just 
P > ImJy in Fig. 1. There appears to be an additional modulation in the whole sample, 
which is less obvious that with P > ImJy, suggesting that this effect is only important for 
the sources with low polarized flux. It appears that there is a modulation with period 180° 
centred around a — 30° which increases the peak at 45° relative to that at —45° which 
is represented by the cos 2 (a — ckq) term. We have illustrated these three effects and their 
combination in Fig. 4 for A = 0.06, B = 0.02, C = 0.02 and ciq = 30°; parameters which 
were decided by trial-and-crror. It is clear that their combination looks very similar to the 
histogram presented in Fig. 1. 

3 SYSTEMATIC EFFECTS ON POLARIZATION MEASUREMENTS 

The model presented in the previous section is purely phenomenological. In this section we 
consider possible real systematic effects which could lead to the biases observed in the data. 
Remarkably, there are reasonable effects which would lead to each of the terms in (2). 

3.1 CLEAN bias 

CLEAN bias is a known phenomenon related to use of the CLEAN algorithm (Hogbom 
1974) on radio synthesis data. Its effect on the data from the NVSS, and an explanation 
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Figure 4. A simple phenomenological model for the observed biases: (bottom left) JVi = 0.06 cos 8a; (top left) JV2 = 0.02 cos 4a; 
(bottom right JV3 = 0.02 cos 2(a - 30°); (top right) the combination of all three N = Ni + N2 + N3. 

for its -working, is discussed in some detail by Condon et al.(1998). Becker et al. (1995) also 
discuss its effect on the FIRST data. 

CLEAN bias produces a reduction in the central flux density of sources by scattering 
flux from the source over a ■wide area. This reduction becomes progressively more severe as 
more CLEAN cycles are used. In the case of the NVSS survey, Condon et al. (1998) state 
that the average clean bias in the NVSS images is approximately —0.3 mJy beam~^. Such 
a bias is relatively benign in the case of total intensity, since it just leads to a small bias in 
the peak flux density, -which is only important for ■weak sources and can be corrected for. 

Its effects on Stokes Q and U maps is much more significant because of (i) the gener- 



NVSS polarization position angles 9 
ally low levels of polarization in radio sources and (ii) Q and U can be both positive and 
negative. The CLEAN bias effect, which is to reduce the magnitude toward zero, can have 
a significant effect on distribution of a even for very small amounts of CLEANing. Since it 
acts proportionately more on either Q or U, which ever is lower. Less information is provided 
on the CLEANing of the Q and U data by the NVSS team; they state that the images are 
"lightly cleaned" since each of the images contain only a few sources (Condon et al. 1998), 
suggesting also that no correction is made for clean bias. 

We model the effects of CLEAN bias by assuming that the observed Stokes parameters 
{Qohs, Uohs) are related to their true values (Qtrue; t^true) by the addition or subtraction 
of ec(> 0) so as to reduce the measured value toward zero. That is, if Qtrue > then 
Qobs = Qtrue - ec, whereas if Qtrue < -ec then Qobs = Qtrue + ec- If -ec < Qtrue < ec then 
Qobs = 0. An equivalent procedure is also used for U. The value of ec can be related to the 
number of times CLEAN subtracts a component from the measured values of Q and U, Nc- 
Using simulations with aQ = au = 300/iJy we have deduced that ec ~ 150/xJy logio(-/Vc); it 
is likely that the coefficient is a function of the noise level. 

In order to investigate the effects of CLEAN bias on the histogram of a, we have produced 
simulated catalogue of sources with similar properties to that of the NVSS. In particular 
we have produced a sample of 500000 sources with polarized fluxes between 0.4mJy and 
lOOmJy drawn from a distribution with dN/dP oc which have polarization angles with 
a uniform distribution. The true values of the Stokes parameters have the effects of CLEAN 
bias included as described above and have Gaussian random noise of 0.3mJybeam~^ added. 
The resulting distribution of position angles is presented in Fig. 5 for ec = 15//Jy and 
300// Jy which have been chosen to illustrate the effect. Both are non- uniform with peaks 
at multiples of 45°, and that with ec = 15//Jy appears to have an amplitude capable of 
explaining the one of the observed effects (oc cos 8a) in the NVSS data. Effects of this kind 
are to be expected since we have already pointed out that CLEAN biases both Q and U 
toward zero, which correspond to multiples of 45° in polarization position angle. 

3.2 Multiplicative and Additive offsets 

Although the bias toward multiples of 45° is the most significant effect in the observed 
data, we have also identified two other possible effects which are oc cos4q; and oc cos 2 (a — 
CKo) respectively. In this section we will show that these effects can be modelled by small 



10 R.A. Battye et al 


























i 



Figure 5. Histogram of a for a simulated NVSS-like survey containing 500000 sources with the effects of CLEAN bias included, 
(left) tQ = 15/xJy and (right) eQ = 300/xJy. In both cases there is a significant bias toward values of a which are multiple of 
45° as is seen in the histogram of the observed data. The histogram can be modelled hy N cos 8a. The errorbars quantify 
the random Poisson error on each bin. 



multiplicative and additive offsets in the data. In particular we will assume that, ignoring 
for the moment CLEAN bias, 



Qobs — (1 + "^Q)Qtrue + Q-Q 



obs 



(1 + mu)Utrue + CLU 



(3) 



that is, the multiplicative offsets are parameterized by rriQ and m^/, and the additive offsets 
by aq and au- 

Errors in polarization calibration are associated with the determination of instrumental 
polarizations for the individual telescopes. The effect of a slight miscalibration typically 
appears as a rotation of points in the (Q, U) plane around a fixed point (Qoj Uq). For a large 
survey with a large number of observing epochs, the effect will be a combination of many 
such offset rotations. For simplicity, we have modelled this complex effect as a combination 
of multiplication and addition in the (Q, U) plane as described above. 

One can compute the systematic error in a due to these two effects to first order 



5a = ^(^[7 ~ itT'q) sin 4a 



2P 



sin 2 (a — ao) 



(4) 



where a is the true polarization angle, P is the polarized flux and tan2ao = clu/O'Q- Convert- 
ing this expression into the corresponding effect on the histogram requires some thought. Let 
us consider which direction the value of a is offset due to the first term which represents the 
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multiplicative offset. There are fixed points at —90°, —45°, 0°, 45°, 90°. If mu—rriQ > then 
—45° and 45° are attractive with the others repulsive. The opposite is true for mu —mq < 0. 
As a general rule we find that each of the trigonometric functions in (4) need to have their 
phase modified by 90°, that is. 



Prom this we see that the multiphcative offset can give rise to an effect oc cos4q; and the 
additive offset can give rise to one oc cos 2 (a — ao). Moreover, we see that the effects due to 
the additive offsets are suppressed by 1/P and hence they would not show up for large P. 

We note that in the case of the multiplicative offset it is really the ratio i? = (1 + 
m[/)/(l + mg) which is relevant and it is this which we will use in the subsequent fitting. 
However, in deriving (4) we have made the assumption that mu and mg are small in which 
case ~ 1 + mu — mq. 

In order to confirm our understanding of these effects we have performed the same sim- 
ulations as discussed in the previous section, but instead of including the effects of CLEAN 
bias, we have introduced a multiplicative offsets of mq = 0.01 and mq = 0.1, while keeping 
mu— 0-0, as well as additive offsets aq — au — 5/xJy and aq — au — 50/xJy. The results of 
these four simulations are presented in Figs. 6 and 7, and these appear to be compatible 
with the prediction (5). 

3.3 Fitting to the observed biases in NVSS 

From the discussion of the preceding sections it appears that a combination of CLEAN bias 
plus multiplicative and additive offsets can lead to histograms qualitatively similar to those 
observed in the NVSS data. Therefore, its seems sensible to establish which values of the 
parameters discussed above best represent the observed biases. We will make an assumption 
in doing this that there are global, flux independent values for these parameters. This may 
not be necessarily true since the amount of CLEANing may vary from field to field, and 
we have not specified the origin of the offsets which could also be different across fields. 
Moreover, we arc presuming that the true values of a arc drawn from a uniform distribution. 

Using the simulation procedure discussed in section 3.1 we have fitted for the parameters 
ec, R = {1+mu) / {1+mq), aq and au in order to give the best fit to the histograms of NVSS 



N = -{mu - mq) sin(4Q; - 90°) 




sin[2(Q; - ao) - 90°] 




(5) 
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Figure 6. Histogram of a for a simulated NVSS survey containing 500000 sources with the effects of multiplicative offsets 
included. (left) niQ = 0.01, mu = 0.0 and (right) mq = 0.1, my = 0.0. The histogram has a bias which can be modelled by 
N oc cos4q. The bias is more visible in the larger case, but the magnitude in the smaller case appears to be similar to that 
required to be added to the CLEAN bias effect to model the observed data. The errorbars quantify the random Poisson error 
on each bin. 
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Figure 7. Histogram of a for a simulated NVSS survey containing 500000 sources with the effects of additive offsets, (left) 
0.Q = 0.U = S/^Jy and (right) aQ = ajj = 50/iJy. The histogram has a bias which can be modelled by TV oc cos2(o — 45°). The 
bias is much more visible in the larger case, but the magnitude of the smaller case appears to be similar to that required to 
be added to the CLEAN bias and multiplicative offsets to explain the observed data in the regime where the polarized flux is 
low. The errorbars quantify the random Poisson error on each bin. 
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Figure 8. Four histograms of the NVSS polarization angles split by polarized flux: (top left) P < IJy; (top right) P = 1 — lOmJy; 
(bottom left) P = 10 — lOOmJy; (bottom right) P > lOOmJy. Included also are a single simulation based on the best fitting 
parameters discussed in the text. 

polarization position angles. We divided the data into separated histograms: P <lmJy, 
1 — lOmJy, 10 — lOOmJy and >100mJy; each with 36 bins of width 5°. The best fits were 
obtained for ec = 14/iJy, R = 0.986, oq = S.lyuJy and au = 6.5/iJy. This value of ec 
corresponds to ~ 1 CLEAN component subtracted per source which seems reasonable at 
least for the weak sources. Note the extremely small offsets applied; these are many times 
smaller than the noise level and could not be detected without such a large sample of objects. 

In general, the fit is good; over all four histograms together, x^/dof(degree of free- 
dom) = 1.60. However, the two higher flux-density histograms are less well fit probably since 
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they contribute proportionately less to the as there are fewer sources in the bins. We 
have attempted to use a flux-dependent value of ec to adjust the level of CLEAN bias and 
hence improve the fit, but with limited success. In practice, the situation is likely to be more 
complicated than our model in this range because of the added effects of self-calibration 
which has been applied during the mapping process of the high flux-density sources. 



3.4 An attempt at correction of the NVSS polarization data 

We have attempted to use the fltted parameters deduced in the previous section to correct 
the NVSS data for the biases in the histogram of a. Such a procedure is fraught with 
danger and should not be taken too seriously since there are two obvious problems: (i) the 
prescription for including CLEAN bias that we have used is not invertible in the sense that 
if IQtruel < then Qohs — 0, and similarly for U ; (ii) only the ratio R can be directly derived 
from the observed data and hence any choice of parameters with ttiq — mu ~ 0.015 will lead 
to the same histogram; (iii) the assumption that the parameters are global and independent 
of polarized flux is unlikely to be completely true. Nonetheless it is an interesting exercise 
since an improvement in the properties would suggest that the model fits source by source, 
rather than just being a pattern spotting exercise on the histograms. 

In order to mitigate against < ec; we have chosen to just perform this "correc- 

tion" on sources with P > ImJy. Moreover, an examination of the two higher flux histograms 
in Fig. 8 suggests that if one also constrains P < 5mJy, one might expect to be at least 
partially successful. We have chosen mg = 1.052 and mu = 1.036 which satisfy the above 
constraint. The resultant histogram is presented in Fig. 9; before the correction the x^/dof 
is ~ 9 and afterwards it is ~ 4 which suggests that, although the procedure we have used is 
obviously not the complete story, there is some reason to believe that it is along the right 
lines. 

We note that although the improvement in the x^/dof is more than a factor of two, 
the change to each of the individual polarization angles is very small. We flnd that the 

((Aa)^)-^/^ pa 0.3° with very few angles changing by more than 2°. It quite interesting to see 
that such small systematic changes in the values of a can lead to such a large effect on the 
xVdof. 
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Figure 9. Histograms of a before (bottom) and after (top) the "correction" in the range Imjy < P < 5mjy for = 14/xJy, 
niQ = 1.052, rrajy = 1.036, aq = 5.1/xJy and ajj = 6.5/imJy. 



4 DISCUSSION AND CONCLUSIONS 

There is a striking systematic in the histogram of a for measurements presented in the 
NVSS catalogue. We have argued that this is as result of a number of very small effects 
that are substantially below the noise level for each source, but which lead to observable 
effects of high significance when "summed" up over the whole catalogue. In attempting to 
make a correction for the effect we find that very small changes in the value a can lead to 
a more uniform histogram. We suggest that CLEAN bias is the chief systematic effect since 
it preferentially selects sources with a = —90°, —45°, 0°, 45° and 90°. Rather than making 
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the correction discussed in section 3.4 which we have aheady pointed out has a number 
of flaws when considered source-by-source, the effects of CLEAN bias can be avoided by 
extracting the sources in the visibihty plane. This was how the final values of Q and U were 
extracted for the JVAS/CLASS survey (Jackson et al 2007) during the analysis of which we 
first became aware of such an effect, albeit for a much smaller number of sources. Adopting a 
visibility plane approach for the NVSS would be much more difficult since there are typically 
a large number of sources in each field, some with complex structure, which was not the 
case in JVAS/CLASS. We emphasize that the effects which we are discussing are generally 
<< 2° for any given source and that using the NVSS catalogue for normal astrophysical 
purposes is unlikely to lead to a substantial error unless the whole catalogue is being used 
statistically to reduce the random errors. 

We note that the kind of biases which we have highlighted may serve as a lesson for the 
future. It appears that before making strong claims on the basis of statistical analysis of 
polarization, data should first be carefully assessed for the kind of effects discussed here. This 
is particularly relevant since new instruments are being designed to perform observations of 
very large numbers of polarized radio sources to study magnetic fields, and of the cosmic 
microwave background (CMB) to search for gravitational waves via the B-mode signatures. 
Both these will require exquisite control of systematics. 

As a final point we discuss how this kind of systematic might affect future high preci- 
sion observations of polarization. Let us first consider a point source survey with the same 
systematic biases (that is, the same values of ec, mq, rriu, aq and au) as in the NVSS, 
but with 2 X 10^ sources between Pmin = 10// Jy and Pmax — lOOmJy with a noise level of 
(yq — — lA*Jy such as might be relevant to the SKA. In such a survey, which is com- 
patible with the distribution of the NVSS observed at higher flux densities, the polarization 
would be detected at SjN > 10 for all sources. It has been shown that this would allow 
one to constrain the rotation measure (RM) to ±5radm~^ using a survey at A = 21cm 
and AA/A — 0.25 (Beck and Gaensler, 2004). A simulated histogram of position angles is 
presented in the left hand column of Fig. 10 showing extremely strong biases, with peaks 
at 0° and 45° which are a factor ~ 3 higher than the rest of the histogram. Such an ef- 
fect would clearly be catastrophic for the extraction of science from such a survey, but is 
probably unrealistic since the levels of CLEAN bias and the additive offsets are likely to be 
defined by the noise level, or possibly the fiux of the source under consideration. We have 
also included in the right hand column of Fig. 10 the results of a simulation with the values 
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Figure 10. Histogram of a for a simulated SKA survey containing 2 X lO'^ sources with systematics included, (left) ec = 15/iJy, 
mg = 0.01, mu = 0.0, aq = 5/iJy and ajj = 6.5/iJy as was found to represent those found in the NVSS survey (right) 
ec = 50njy, mq = 0.01, mjj = 0.0, ag = 17njy and au = 22njy which correspond to a scaling of the noise level expected for 
the SKA survey, relative to the NVSS survey. 

of ec, cJq and au scaled down by a factor 300 which is the ratio of the noise level in the 
proposed SKA survey to that in the NVSS. The amplitude of bias in this case is of similar 
order of magnitude to that found in the NVSS although the structure of the peaks in the 
histogram is much better defined due to the increased number of sources in each bin. 

It is interesting to speculate on how the smaller of these two biases might affect the 
extraction of science from such a survey. If the bias is similar in character to that in the 
NVSS then ((Aa)^) ^ 0.3°, that is, there is a offset, dependent on the intrinsic position 
angle with this RMS. This will lead to an offset in the measured RM and its error. The noise 
error on each position angle will be around 5° corresponding to the the S/N > 10 which 
is much larger than the off set /increase in the noise. However, it is planned for the sets of 
rotation measures to be used together to reduce the random errors and deduce information 
about magnetic fields. If ~ 300 sources were added together then the resultant random errors 
would be comparable to the offset. 

The other area where high precision measurements of polarization are planned is the 
CMB. There position angle accuracies << 0.3° will be required to probe primodial gravi- 
tational waves corresponding to scalar-to-tensor ratios, r ~ 10^'^ which might be possible 
within planned polarization satellites. In this application it is unlikely that CLEAN bias will 
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be a problem since power spectrum estimates will be made directly from the "dirty" images 
and hence do not use the CLEAN algorithm, but the multiplicative and additive offsets are 
likely to be an issue. We plan to investigate whether the kind of techniques used here can 
be used to search for subtle systematic effects in the measured polarization. 
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